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Abstract — Cognitive  radio  (CR)  is  an  enabling  technology  for 
efficient  use  of  available  spectrum  and  promises  unprecedented 
flexibility  in  multi-hop  wireless  networking.  This  paper  explores 
networking  related  issues  associated  with  CRs.  Specifically,  we 
consider  how  to  maximize  the  rates  of  a  set  of  user  communi¬ 
cation  sessions  in  a  multi-hop  CR-based  wireless  network.  Due 
to  potential  interference  at  the  physical  layer,  we  find  that  it  is 
essential  to  follow  a  cross-layer  approach,  with  joint  optimization 
at  physical  (power  control),  fink  (frequency  band  scheduling), 
and  network  (flow  routing)  layers.  We  give  a  mathematical 
characterization  of  this  cross-layer  optimization  problem.  We 
develop  a  centralized  solution  procedure  based  on  the  branch- 
and-bound  framework.  Using  numerical  results,  we  demonstrate 
the  efficacy  of  the  solution  procedure  and  offer  quantitative 
understanding  on  the  joint  optimization  at  different  layers. 

Index  Terms — Cognitive  radio,  multi-hop  wireless  network, 
interference,  cross-layer  optimization. 

I.  Introduction 

Cognitive  radio  (CR)  is  a  revolution  in  radio  technology 
that  is  enabled  by  advances  in  RF  design,  signal  processing, 
and  communications  software  [14],  It  promises  unprecedented 
flexibility  in  radio  communications  and  efficient  use  of  spec¬ 
trum.  The  potential  of  CR  has  been  recognized  by  the  commer¬ 
cial  sector  as  well  as  the  military  (e.g.,  JTRS  program  [10]) 
and  public  safety  communications  (e.g.,  SAFECOM  [15]). 

Our  goal  in  this  paper  is  to  optimize  network  level  perfor¬ 
mance  of  multi-hop  CR  networks.  It  is  now  well  understood 
that  network  performance  for  such  networks  is  tightly  coupled 
with  lower  layer  behaviors  [17].  For  instance,  maximizing  user 
throughput  at  the  network  level  not  only  depends  on  flow 
routing,  but  also  depends  on  the  algorithms  at  link  layer  (e.g., 
frequency  band  assignment)  and  physical  layer  (e.g.,  power 
control).  As  a  result,  an  optimal  solution  at  the  network  level 
must  be  developed  with  joint  consideration  of  multiple  layers. 

Recently,  there  is  a  growing  interest  on  gaining  under¬ 
standing  on  multi-hop  CR  networks  (see,  e.g.,  [17],  [18]). 
However,  most  of  these  work  are  based  on  the  so-called 
“protocol  interference  model”  [8].  Under  such  model,  the 
notions  of  transmission  range  and  interference  range  are  used 
to  determine  the  feasibility  of  successful  transmission  and  the 
existence  of  interference.  It  is  now  well  understood  that  such 
binary  decision  on  interference  modeling  has  its  limitation. 
On  the  other  hand,  the  so-called  “physical  model”  is  widely 
accepted  as  an  accurate  characterization  of  interference.  Under 
physical  model,  a  transmission  is  successful  if  and  only  if 
signal-to-interference-and-noise-ratio  (SINR)  at  the  intended 


receiver  exceeds  a  certain  threshold  so  that  the  transmitted 
signal  can  be  decoded  with  an  acceptable  bit  error  rate 
(BER).  Further,  capacity  calculation  is  based  on  SINR  (via 
Shannon’s  formula),  which  takes  into  account  of  interference 
due  to  simultaneous  transmission  at  other  nodes.  Unfortu¬ 
nately,  although  physical  model  is  accurate,  there  is  much 
difficulty  in  carrying  out  analysis  with  physical  model  due 
to  the  computational  complexity  it  involves,  particularly  when 
it  comes  to  cross-layer  optimization  in  a  multi-hop  network 
environment. 

In  this  paper,  we  investigate  networking  problem  for  multi¬ 
hop  CR  networks.  We  employ  the  physical  model  for  in¬ 
terference  modeling  and  study  the  problem  via  cross-layer 
optimization  approach.  In  particular,  we  consider  how  to 
maximize  the  rates  of  a  set  of  user  communication  sessions, 
with  joint  consideration  at  physical  layer  (via  power  control), 
link  layer  (via  frequency  band  scheduling),  and  network  layer 
(flow  routing).  We  give  a  mathematical  characterization  for 
these  layers  and  formulate  the  problem  into  a  mixed  integer 
nonlinear  program  (MINLP).  We  develop  a  centralized  solu¬ 
tion  to  this  complex  optimization  problem  based  on  branch- 
and-bound  (BB)  framework  and  a  reformulation-linearization 
technique  (RLT)  [16].  The  solution  we  develop  is  guaranteed 
to  be  within  a  factor  of  (1  —  e)  from  the  optimum,  where  e  is 
a  small  parameter  reflecting  our  desired  accuracy. 

The  remainder  of  this  paper  is  organized  as  follows.  In 
Section  II,  we  give  a  mathematical  characterization  of  power 
control,  scheduling,  and  routing  for  multi-hop  CR  networks. 
We  also  present  the  problem  formulation  in  this  study.  In 
Section  III,  we  present  a  solution  based  on  branch-and-bound 
framework.  Section  IV  presents  numerical  results  for  the  cross¬ 
layer  solution.  Section  V  reviews  related  work  and  Section  VI 
concludes  this  paper. 

II.  Mathematical  Models 

We  consider  a  CR-based  ad  hoc  network  with  a  set  of  nodes 
A f.  For  a  node  i  £  AT,  the  set  of  available  frequency  bands 
A/1,  depends  on  its  location  and  may  not  be  identical  to  the 
available  frequency  bands  at  other  nodes.  We  assume  that  the 
bandwidth  of  each  frequency  band  (channel)  is  W.  Denote 
Ai  the  set  of  all  frequency  bands  present  in  the  network,  i.e., 
M  =  IW  A ii.  Denote  M-ij  =  which  is  the 

set  of  frequency  bands  that  is  common  on  both  nodes  i  and 
j  and  thus  can  be  used  for  transmission  between  these  two 
nodes.  In  the  rest  of  this  section,  we  present  mathematical 
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TABLE  I 
Notation. 


Symbol 

Definition 

N 

The  set  of  nodes  in  the  network 

Mi 

The  set  of  available  bands  at  node  i  E  J\f 

M 

The  set  of  frequency  bands  in  the  network 

Mij 

The  set  of  frequency  bands  on  link  i  — >  j 

W 

Bandwidth  of  a  frequency  band 

C 

The  set  of  active  user  communication  sessions 

s(l),  d(l) 

Source  and  destination  nodes  of  session  Z  E  C 

r(l) 

Minimum  rate  requirement  of  session  l 

K 

Rate  scaling  factor  for  all  sessions 

Pmax 

The  maximum  transmission  power  at  a  transmitter 

V 

Ambient  Gaussian  noise  density 

9ij 

Propagation  gain  from  node  i  to  node  j 

a 

The  minimum  required  SINR 

rj~m 

1  i 

The  set  of  nodes  that  can  transmit  to  (and  receive 
from)  node  i  on  band  m 

T, 

The  set  of  nodes  that  can  transmit  to  (and  receive 
from)  node  i 

'Tm 

The  set  of  nodes  that  may  make  interference  on 
band  m  at  node  j 

Xij 

Binary  indicator  to  mark  whether  or  not  band  m 
is  used  by  link  i  — >  j 

fail ) 

Data  rate  for  session  l  on  link  i  — >  j 

Q 

The  number  of  transmission  power  levels 

nVt 

The  transmission  power  level  from  node  i  to  node 
j  on  band  m 

j.m 

The  transmission  power  level  at  node  i  on  band  m 

Qm 

_ _ 

The  SINR  from  node  i  to  node  j  on  band  m 

Notation  for  branch-and-bound  procedure 

£ 

A  small  positive  constant  reflecting  our  desired 
accuracy  in  the  final  solution 

nz 

The  set  of  all  possible  solutions  in  problem  z 

LBZ,  UB, 

The  lower  and  upper  bounds  of  problem  z 

Ipz 

The  local  search  solution  for  problem  z 

LB,  UB 

The  maximum  lower  and  upper  bounds  among  all 
problems 

i’e 

The  (1  —  e)  optimal  solution 

characterization  of  each  layer  in  a  multi-hop  CR  network. 
Table  I  lists  all  notation  in  this  paper. 


A.  Scheduling  and  Power  Control 

Scheduling  for  transmission  at  each  node  in  the  network 
can  be  done  either  in  time  domain  or  frequency  domain.  In 
this  paper,  we  consider  scheduling  in  frequency  domain  in  the 
form  of  assigning  frequency  bands  (channels).  To  maximize 
the  capacity,  there  may  still  be  concurrent  transmissions  within 
the  same  channel  (and  thus  interference). 

Denote 


{1  If  node  i  transmits  data  to  node  j  on  band  m, 
0  otherwise. 


(1) 


Due  to  interference,  a  node  i  can  use  a  band  m  for  transmitting 
to  one  node  j  or  receiving  from  one  node  k.  That  is. 


E  x%  +  E  xv  ^ 1 .  (2) 

feer™  jeTx1 


where  7Jm  is  the  set  of  all  possible  nodes  that  node  i  can 
transmit  to  (and  receive  from)  on  band  m  in  the  network. 

For  power  control,  we  assume  that  the  transmission  power 
at  a  node  can  be  tuned  to  a  finite  number  of  levels  between  0 
and  Pmax.  To  model  this  discrete  version  of  power  control,  we 
introduce  an  integer  parameter  Q  that  represents  the  total  num¬ 
ber  of  power  levels  to  which  a  transmitter  can  be  adjusted,  i.e.. 


0,  ^Pmax,  §Tjnax,  Pmax-  Denote  q™  e  {0, 1, 2,  •  •  • ,  Q}  the 
integer  power  level.  Clearly,  when  node  i  does  not  transmit 
data  to  node  j  on  band  to,  q™  should  be  0.  Under  the 
maximum  allowed  transmission  power  level  Q,  we  have 


<  Q 
=  0 


If  x%  =  1, 

otherwise. 


With  joint  consideration  of  x™  and  q™,  the  above  relationship 
can  be  re-written  as 


q?,  <  Qxi 


(3) 


which  shows  that  the  coupling  relationship  between  power 
control  and  scheduling. 

As  discussed  earlier,  to  maximize  capacity,  there  may  be 
concurrent  transmissions  by  different  nodes  on  the  same  band. 
Under  physical  model,  a  transmission  is  successful  if  and  only 
if  the  SINR  at  the  receiving  node  exceeds  a  certain  threshold, 
say  a.  Mathematically,  for  a  transmission  from  node  i  to 
node  j  on  band  to,  when  there  is  interference  from  concurrent 
transmissions  on  the  same  band,  the  SINR  is 


Si  A  = 


n  Sii-P 


vW  +  E 


fceAf  ELr, 
9ij9 Tj 


9kj-$ 


Pm 


rjWQ 

Pmax 


Z^keAf  Z^ 


h£T™  9kj9kh 


where  tj  is  the  ambient  Gaussian  noise  density,  (jtJ  is  the 
propagation  gain  from  node  i  to  node  j. 

Based  on  our  model  for  scheduling,  if  x ^  =  1,  then  x'ki  =  0 
for  k  £  T™  and  x™-  =  0  for  k  £  7”m.  As  a  result,  by  (3), 
=  0  and  q™-  =  0  Then  we  have 


s  a  = 


rjWQ 

/'max 


9,.,'ir, 

2-ikeN  A-iheT” 


9kj9kh 


Denote  C,  We  have 

9ijQij 


Si  A  = 


_  ,n 

Pmax  ;T«  9kjlik 

Note  that  this  SINR  computation  also  holds  when  <■/[' 
i.e.,  when  there  is  no  transmission  from  node  i  to  node  j  on 
band  to. 

Recall  that  under  physical  model,  a  transmission  from  node 
i  to  node  j  on  band  m  is  successful  if  and  only  if  s?j  > 
Then  by  (1),  we  can  couple  x™  and  s"i-  as  follows. 


rjWQ 


(4) 

=  0, 


If  s%  >  a. 


)  otherwise. 

which  can  be  written  into  the  following  equivalent  relationship. 


sTj  >  ax™ 


B.  Routing 

In  a  multi-hop  ad  hoc  network,  we  assume  there  is  a  set  of 
£  active  user  communication  (unicast)  sessions.  Denote  s(l) 
and  d(l)  the  source  and  destination  nodes  of  session  l  £  £ 
and  r(l)  the  minimum  rate  requirement  (in  b/s)  of  session  l. 
In  our  study,  we  aim  to  maximize  a  scaling  factor  K  for  all 


session  rates.  That  is,  what  is  the  maximum  factor  K  such 
that  a  rate  of  K  ■  r(l)  can  be  transmitted  from  s(l)  to  d(l)  for 
each  session  l  £  C  in  the  network. 

To  route  these  data  flows  from  its  source  node  to  destination 
node,  multi-hop  relaying  is  necessary,  due  to  limited  transmis¬ 
sion  power  at  each  node.  Further,  for  optimality  and  flexibility, 
it  is  desirable  to  allow  flow  splitting  and  multi-path  routing. 
This  is  because  a  single  path  flow  routing  for  a  session  is 
overly  restrictive  and  is  unlikely  to  guarantee  optimal  solution. 

Mathematically,  this  can  be  modeled  as  follows.  Denote 
fij(l)  the  data  rate  on  link  (i.  j)  that  is  attributed  to  session 
l,  where  i  £  Af.j  £%  =  UmeAt  ^m-  If  node  i  is  the  source 
node  of  session  l,  i.e.,  i  =  s(l),  then 

E  fad)  =  r(l)K  .  (5) 

If  node  i  is  an  intermediate  relay  node  for  session  l,  i.e.,  i  ^ 
s(l)  and  i  ^  d(l ),  then 

j^s{l)  k^d{l) 

ei  =  ei  fki (o  •  (6) 

jeT  keTi 

If  node  i  is  the  destination  node  of  session  l,  i.e.,  i  =  d(l), 
then 

E  /w(0  =  r(l)K  ■  (7) 

k^Ti 

It  can  be  easily  verified  that  once  (5)  and  (6)  are  satisfied,  (7) 
must  also  be  satisfied.  As  a  result,  it  is  sufficient  to  list  only 
(5)  and  (6)  in  the  formulation. 

In  addition  to  the  above  flow  balance  equations  at  each  node 
i  £  J\f  for  session  l  £  C ,  the  aggregated  flow  rates  on  each 
radio  link  cannot  exceed  this  link’s  capacity.  For  a  link  i  — >  j, 
we  have 

E  /«(0<  E  Wlog2(l  +  s%).  (8) 

1£jC  tyi  £  J^A.  i  j 

The  constraint  in  (8)  further  illustrates  the  coupling  relation¬ 
ship  among  flow  routing,  power  control,  and  scheduling. 

C.  Problem  Formulation 

Putting  together  all  the  constraints  for  scheduling,  power 
control,  and  flow  routing,  we  have  the  following  complete 
problem  formulation. 

Max  K 

s.t.  E  x™i  +  E  —  I  ( i£jV,m£Mi ) 

q?j~Qx?j<Q  (i£AT,m£Mi,j£  %m )  (9) 

E  Qij  ~  tT  =  0  (i  £  JV,m  £  Mi ) 

WD 

Tj VV  (q/  m  \  A  ,m  m  m 

P  T  /  /  fjk  'i  /'  k  $ij  Qij  Qij  u 

( i£U,m£  Mi,j  £  T"1) 

(i  £  Af,m  £  Mi,j  £  %m)  (10) 


''I'  ^  f\ 

-  si:i  <  0 


s<J)^j,d(t)^i 

E  Ai(0-  E  w\og2(i  +  s™)<0 

l££ 

(■ i  €  A r,j  £  %) 

E  fijil)  ~  r{l)K  =  0  (l  £  C,i  =i  s(l)) 

jeT 

j^s(l)  k^d(l) 

E  E  (l£C,i£tf,i?s(l),d(l)) 

jeT  ierk 

x?j  £  {0,1  },qTj  £  {0, 1,  2,  •  •  • ,  Q},  tT ,  s™  >  0 

(i  £  N ,  m  £  Mi,j  £  T/n) 

K,  fijil)  >0  {l  £  C,i  £  N ,i  =£  d(l),j  £%,j  ^  s(l ))  , 

where  Q,r],W,a,  Pmax,  gij,  and  r(l)  are  all  constants  and 
K,  x™ ,  q™,  t™,  s™,  and  fij(l)  are  all  optimization  variables. 
This  formulation  is  a  mixed  integer  non-linear  program 
(MINLP),  which  is  NP-hard  in  general  [6], 

III.  A  Centralized  Solution 

A.  Oven’iew 

For  the  complex  MINLP  problem,  we  employ  the  so-called 
branch-and-bound  framework  [13]  to  develop  a  solution.  Un¬ 
der  branch-and-bound,  we  aim  to  provide  a  (1  —  e)-optimal 
solution,  where  e  is  a  small  positive  constant  reflecting  our 
desired  accuracy  in  the  final  solution. 

To  start  with,  branch-and-bound  employs  some  relaxation 
technique  to  obtains  a  linear  relaxation  for  the  original  prob¬ 
lem.  The  solution  to  this  relaxed  problem  provides  an  upper 
bound  (UB)  to  our  objective  function.  With  the  relaxation 
solution  as  a  starting  point,  branch-and-bound  uses  a  local 
search  algorithm  to  find  a  feasible  solution  to  the  original 
problem,  which  provides  a  lower  bound  (LB)  for  the  objective 
function.  If  the  obtained  lower  and  upper  bounds  are  close  to 
each  other,  i.e.,  LB  >  (1  —  e)UB,  then  the  current  feasible 
solution  is  (1  —  e)  optimal  and  we  are  done. 

Otherwise,  branch-and-bound  replaces  the  original  prob¬ 
lem  with  two  sub-problems,  say  problem  1  and  problem  2. 
This  is  accomplished  by  choosing  an  appropriate  variable  for 
partitioning  and  dividing  the  range  of  the  partition  variable 
into  smaller  ranges.  For  problems  1  and  2,  branch-and-bound 
performs  relaxation  and  local  search  on  each  problem  and 
we  obtain  upper  bounds  UB i  and  U IU  and  lower  bounds 
LB i  and  LB2  for  problems  1  and  2,  respectively.  Since  the 
relaxations  in  problems  1  and  2  are  both  tighter  than  that  in 
the  original  problem,  we  have  max{UBi,  UB2}  <UB  and 
max{LB\,  LB2}  >  LB.  Then  the  upper  bound  of  the  original 
problem  is  updated  as  UB  =  max{UBi,  UB2}  and  the  lower 
bound  is  updated  as  LB  =  max{LBi,  LB2}.  As  a  result, 
we  now  have  smaller  gap  between  UB  and  LB.  Then  we 
either  have  a  (1  — e)  optimal  solution  (if  LB  >  (l  —  e)UB)  or 
choose  a  problem  with  the  maximum  upper  bound  and  perform 
partition  for  this  problem. 

An  important  property  of  branch-and-bound  is  that  we  may 
remove  some  problems  from  further  consideration  during  the 
iterations.  In  particular,  if  we  find  a  problem  z  with  LB  > 
(1  —  s)UBz,  then  we  can  conclude  that  this  problem  cannot 


Fig.  1.  A  convex  hull  for  y  =  \nx. 

provide  much  improvement  on  LB  and  we  can  thus  remove 
this  problem  from  further  consideration. 

In  the  rest  of  this  section,  we  present  the  key  components  in 
the  branch-and-bound  framework,  which  are  problem  specific 
and  far  from  trivial. 


Local  Search  Algorithm 

Initialization: 

1.  Set  q =  {q™)L  and  x^j  as  0  or  1  if  its  value  set 
only  have  one  element  0  or  1,  respectively. 

2.  Compute  the  requirement  fij  (/). 

Iteration: 

3.  Compute  the  capacity  and  the  ratio  kij  between  the  capacity 
and  the  requirement  for  each  link  i  — *■  j. 

4.  Suppose  link  i  — >  j  has  the  smallest  kij.  We  will  try  to 
increase  its  capacity  as  follows. 

5.  If  we  can  increase  q Yj  on  a  used  band  { 

6.  Suppose  band  m  has  the  largest  q ™  among  these  bands. 

7.  Increase  q™j  under  the  constraints  of  q™  <  (q™ )  u  and 

new  kij  <  1.  } 

8.  else,  if  we  can  increase  q™  on  an  available  but  currently 
unused  band  { 

9.  Suppose  band  m  has  the  largest  q™  among  these  bands. 

10.  Increase  q ™  under  the  constraints  of  q <  ( Qij)u  and 

new  kij  <  1. 

1 1.  We  also  need  to  set  x1^  =  1,  xYh  —  0  for  h^Ti,  h^j, 

=  0  for  k  £  Ti.  } 

12.  else  the  iteration  terminates  since  the  smallest  kij  cannot  be 
increased. 


Fig.  2.  Pseudocode  of  proposed  local  search  algorithm. 


constraints. 


B.  Linear  Relaxation 


During  each  iteration  of  the  branch-and-bound  procedure, 
we  need  a  linear  relaxation  technique  to  obtain  an  upper  bound 
of  the  objective  function. 

For  the  polynomial  term  in  the  problem  formu¬ 

lation,  we  apply  a  novel  method  based  on  Reformulation- 
Linearization  Technique  (RLT)  [16].  That  is,  we  introduce 
a  new  variable  replace  by  and  add  RLT 

constraints  on  these  variables.  Suppose  trk  and  s™  are  bounded 
by  {t?)L  <  t?  <  (f T)u  and  (s%)L  <  s™  <  {8%)Ut 
respectively.  Thus,  we  have  [t™  -  (f™)z,]  •  [s™  -  (s[")l]  >  0, 
K-^lH^u-sV*)  >  0,  [(tT)u-tT]-K-(s^)L]  > 


0,  and  [(tZ)u  -  q 

relationships  and  si 
following  RLT  constraints  for  u-lk. 


[(si7)tt  —  sJ7l  >  0.  From  the  above 


iji  — 

relationships  and  substituting  u™k  =  t™s™,  we  have  the 


*. r 


(tf)L 

(t™)u 

(wh 

(t™)u 


„m  j  ( c>m\  xm  „,m  <<  ( xm\  ( „m\ 

sij  +  {sij)L  ’  tk  —  Ujjfc  <  [tk  )L  '  \Sij)L  , 

nm  i  ( Cim\  xm  „,m  \  / j.m\  ( „7n\ 

& ij  T  '  tk  U ijk  —  [t-k  )u  •  \S'ij  )  L  7 

l  (  j.712  „ , m  \  (4 .m\  (  ttm\ 

sij  +  Vsij)u  ■  tk  -  uijk  >(tk)L-  (Sij)u  , 
sij  +  ( sij)u  ■  tk  ~  Uijk  <  (tk  )u  ■  . 


For  the  log  term,  we  propose  to  employ  three  tangential 
supports,  which  is  a  convex  hull  linear  relaxation.  We  first 
analyze  the  bounds  for  1  +  sf!j.  Then,  we  introduce  a  variable 
c™  =  ln(l  +  8%)  and  consider  how  to  get  a  linear  relaxation 
for  y  =  In  x  over  xl  <  x  <  xjj-  This  function  can  be  bounded 
by  four  segments  (or  a  convex  hull),  where  segments  I,  II,  and 
III  are  tangential  supports  and  segment  IV  is  the  chord  (see 
Fig.  1).  In  particular,  three  tangent  segments  are  at  (xl,  lnx^), 
(/?, ln/3),  and  (xjj ,\nxu),  where  [3  =  XL'Xu ln x L )  js 
the  horizontal  location  for  the  point  intersects  extended  tangent 
segments  I  and  III;  segment  IV  is  the  segment  that  joins  points 
(xi,  ln  xl)  and  {xjj,  lnxy).  The  convex  region  defined  by  the 
four  segments  can  be  described  by  the  following  four  linear 


xl  ■  y  -  x  <  XL(lnXL  -  1)  , 

P-y~x<f3{ ln/3— 1)  , 

Xu  ■  y  -  X  <  Xr/(ln  Xu  -  1)  , 

(xu~ XL)y+(hxxL—  ln  xu)x>xu  lnxL  —  Xl  In  Xu  ■ 

As  a  result,  the  non-polynomial  (log)  term  can  also  be  relaxed 
into  linear  constraints. 

Based  on  the  above  linear  relaxation  techniques,  we  can 
relax  the  original  problem  into  a  linear  program  (LP). 


C.  Local  Search  of  Feasible  Solution 

A  linear  relaxation  for  a  problem  z  can  be  solved  in 
polynomial  time.  Denote  the  relaxation  solution  as  if)z,  which 
provides  an  upper  bound  to  problem  z  but  may  not  be  feasible. 
We  now  show  how  to  obtain  a  feasible  solution  if)z  based  on 

i’z- 

Denote  x  and  q  as  the  vector  for  variables  x™  and  qf, 
respectively.  To  obtain  a  feasible  solution,  we  need  to  deter¬ 
mine  the  integer  values  for  x  and  q  in  solution  fz  such  that 
(2),  (4),  (9),  (10)  hold.  All  other  variables  are  based  on  x  and 
q.  Initially,  each  q™  is  set  to  the  smallest  value  (qf-)L  in  its 
value  set  and  x™  is  fixed  as  0  or  1  if  its  value  set  only  has  one 
element  0  or  1,  respectively.  Based  on  these  q™’ s,  we  can  com¬ 
pute  the  capacity  WloS2  (  1  +  riWQ  .  \  -  A:  /  , . ;;  ,  ?11  j 

3  \  i=max  9kjtk  J 

for  each  link  i  — >  j.  The  requirement  on  a  link  i  — >  j 
is  ^ (if  Thus,  we  can  compute  kij,  the  ratio 
between  the  capacity  and  the  requirement.  The  objective  value 
for  the  current  x  and  q  is  I\  ■  min  {kij  :  i  £  A f,j  £  7)}. 
Thus,  we  aim  to  increase  the  minimum  kij.  We  always  try 
to  increase  the  smallest  kij  by  increasing  some  qfj  under  the 
constraint  r/"'  <  ( qf!j)u-  When  we  cannot  further  increase  the 
smallest  kij,  we  are  done.  The  pseudocode  of  this  local  search 
algorithm  is  given  in  Fig.  2. 


TABLE  II 

Location  and  available  frequency  bands  at  each  node  for  a  20-node  5-session  network. 


Node 

Location 

Available  Bands 

Node 

Location 

Available  Bands 

Node 

Location 

Available  Bands 

1 

(0.1, 9.9) 

1,  2,  3,  4,  7,  8,  9,  10 

8 

(22.6,40.9) 

1,  2,  3,  5,  7,  9,  10 

15 

(44.7,  24) 

1,  2,  3,  4,  5,  6,  7,  8,  9,  10 

2 

(29.2,31.7) 

1,  2,  3,  4,  5,  7,  8,  10 

9 

(35.3,  10.3) 

2,  9 

16 

(47.9,43.8) 

1.  3 

3 

(3,31.1) 

1,  4,  5,  6 

10 

(31.9, 19.6) 

1,  2,  3,  4,  5,  6,  7,  8,  9,  10 

17 

(46.4,  16.8) 

1,  7,  9 

4 

(11.8,40.1) 

1,  2,  3,  4,  6,  9,  10 

11 

(28.1,25.6) 

1,  2,  3,  4,  5,  6,  7,  8,  9,  10 

18 

(11.5, 12.2) 

2,  5,  6,  10 

5 

(15.8,9.7) 

1,  2,  3,  5,  6.  8,  9 

12 

(32.3,38) 

1,  8,  9,  10 

19 

(28.2,  14.8) 

4,  5,  6,  7,  8,  9,  10 

6 

(16.3,  19.5) 

3.  5,  6.  8,  9 

13 

(47.2,2.6) 

3,  5,  10 

20 

(2.5, 14.5) 

1,  7,  10 

7 

(0.6,27.4) 

1,  4,  8,  9,  10 

14 

(44.7, 15) 

2,  3,  6,  7,  8 

TABLE  III 

Source  node,  destination  node,  and  minimum  rate  requirement 

OF  EACH  SESSION  IN  THE  20-NODE  5-SESSION  NETWORK. 


Session 

l 

Source  Node 
s(l) 

Dest.  Node 
d(i) 

Min.  Rate  Req. 

r(l) 

1 

16 

10 

9 

2 

18 

3 

1 

3 

12 

11 

4 

4 

13 

17 

3 

5 

15 

6 

2 

D.  Selection  of  Partition  Variables 

Based  on  the  impact  on  the  objective  value,  variables  in  x 
are  more  important  than  variables  in  q.  Thus,  we  should  first 
select  one  of  x  variables  as  the  branch  variable.  In  particular, 
for  the  relaxation  solution  ipz,  the  relaxation  error  of  a  discrete 
variable  x”J  is  minjxJ'Tl  —  x"’},  where  x",-'  is  the  value 

LJ  V.  LJ  ^  L  J  J  LJ 

of  variable  x^j  in  solution  %j)z.  We  choose  an  x™  with  the 
maximum  relaxation  error  among  all  x  variables  and  let  its 
value  set  in  problems  Z\  and  z2  be  {0}  and  {1},  respectively. 
Since  the  value  set  for  this  x("  only  has  one  element,  this  xf] 
can  be  replaced  by  a  constant  in  the  new  problem.  As  a  result, 
some  constraints  may  also  be  removed. 

It  should  be  note  that  we  may  pose  more  limitations  on 
other  variables  based  on  the  new  value  set  of  x”j.  That  is,  if 
the  new  value  set  of  xj*j  is  {0},  then  we  have  qfP  =  0  based  on 
(9).  If  the  new  value  set  of  x™  is  {1},  then  we  have  xff  =  0 
for  h  £  Ti,h  ^  j  and  xj™  =  0  for  k  €  based  on  (2). 

When  none  of  the  x  variables  can  be  partitioned  (i.e., 
each  of  their  value  sets  has  only  one  element),  we  select 
one  of  q  variables  for  partitioning.  In  particular,  in  the  re¬ 
laxation  solution  ipz,  the  relaxation  error  of  g”1  is  min {q™  — 
Ig^J ,  [g™J  +  1  —  g,™},  where  g™  is  the  value  of  variable 
q™  in  solution  ipz.  Assuming  the  value  set  of  g™  in  problem 
z  is  {go,  gi,  •  •  •  ,Qk},  its  value  set  in  problems  Z\  and  z2 
will  be  {g0,  gi,  •  •  • ,  [q%\ }  and  { |_g^J  +1,  [q%\  +  2,  •  •  • ,  qK}, 
respectively.  Again,  based  on  the  new  value  set  of  gj™,  we  may 
impose  additional  limitations  on  other  variables.  In  particular, 
if  the  new  value  set  of  g™  is  {0},  then  we  have  x”j  =  0  based 
on  (10).  If  the  new  value  set  of  g™  does  not  include  0,  then 
we  have  x”j  =  1  based  on  (9). 

Note  that  when  all  possible  partition  variables  in  x  and  q 
can  no  longer  be  partitioned  (i.e.,  all  values  are  assigned),  the 
other  variables  can  be  solved  via  an  LP. 

IV.  Numerical  Results 

In  this  section,  we  present  numerical  results  on  the  proposed 
solution.  Our  goals  are  to  demonstrate  the  efficacy  of  the 


solution  procedure  and  offer  quantitative  understanding  on  the 
joint  optimization  at  different  layers. 

A.  Simulation  Setting 

For  the  ease  of  exposition,  we  normalize  all  units  for 
distance,  bandwidth,  rate,  and  power  based  on  (8)  with  ap¬ 
propriate  dimensions.  We  consider  a  20-node  CR  networks 
with  each  node  located  in  a  50x50  area.  We  assume  there  are 
\AA\  =  10  frequency  bands  in  the  network  and  each  band  has 
a  bandwidth  of  W  =  50.  At  each  CR  node,  only  a  subset 
of  these  bands  is  available.  Table  II  gives  the  details  of  the 
location  of  each  node  and  the  set  of  available  bands  at  each 
node.  We  assume  there  are  5  user  communication  sessions, 
each  with  a  minimum  rate  requirement  within  [1,10].  The 
source  node,  destination  node,  and  minimum  rate  requirement 
of  each  session  are  given  in  Table  III. 

We  assume  the  propagation  gain  is  j,:/  =  cZ~-  1  and  the  SINR 
threshold  a  =  3  [7],  The  maximum  transmission  power  at  each 
node  is  Pmax  =  4.8  •  105  •  rj  •  W.  We  assume  that  power  control 
can  be  done  in  Q  =  10  levels. 

For  our  proposed  branch-and-bound  solution  procedure,  we 
set  e  to  0.1,  which  guarantees  that  the  solution  is  within  90% 
optimal. 

B.  Results  and  Observations 

For  the  20-node  network  with  5  sessions,  the  transmission 
power  levels  on  their  respective  frequency  bands  in  the  final 
solution  are: 

Band  1:  g{3  =  1,  g{6  12  =  7; 

Band  2:  gf  2  =  2; 

Band  3:  g4314  =  2; 

Band  4:  qf  7  =  7,  g|10  =  2; 

Band  5:  q\l  w  =  1; 

Band  6:  qf5  19  =  9; 

Band  7:  g[417  =  1,  g|01  =  1; 

Band  8:  g42  n  =  3; 

Band  9:  qf2  8  =  1,  g?9  6  =  3; 

Band  10:  g{g  20  =  1. 

Note  that  the  same  frequency  band  may  be  used  by  concurrent 
transmissions,  e.g.,  both  node  7  — >  3  and  node  16  — >  12  are 
transmitting  on  band  1.  To  minimize  interference,  our  solution 
has  placed  these  concurrent  transmissions  sufficiently  apart 
and  set  the  optimal  transmission  power  less  than  the  maximum. 

Figure  3  shows  the  routing  topology  in  the  final  solution. 
The  flow  rates  are 

Session  1:  /2,10(1)  =  103.30,  /8>2(1)  =  103.30,  /n,i0(l)  = 
15.86,  /12,8(  1)  =  103.30,  /i2,ii(l)  =  15.86,  /i6,i2(l)  = 


Fig.  3.  The  routing  topology  for  the  20-node  5-session  network. 

119.16; 

Session  2:  /1)7( 2)  =  13.24,  /7>3( 2)  =  13.24,  /i8)20(2)  = 
13.24,  /20,i(2)  =  13.24; 

Session  3:  /i2,n(3)  =  52.96; 

Session  4:  /i3’i4(4)  =  39.72,  /i4.i7(4)  =  39.72; 

Session  5:  /15’19(5)  =  26.48,  /i9,6(5)  =  26.48. 

We  can  see  that,  to  maximized  the  achieved  capacity,  multi- 
path  routing  are  used  for  session  1. 

Under  this  solution,  the  achieved  data  rate  for  sessions  1 
to  5  is  119.16,13.24,52.96,39.72,26.48,  respectively,  which 
corresponds  to  a  scaling  factor  of  13.24. 

V.  Related  Work 

There  has  been  substantial  research  efforts  on  multi-hop 
wireless  networks  based  on  the  protocol  interference  model 
(see,  e.g.,  [1],  [9],  [11]).  The  controversy  surrounding  (or  argu¬ 
ments  against)  the  protocol  model  is  that  a  binary  decision  of 
whether  or  not  interference  exists  (based  on  interference  range) 
does  not  accurately  captures  physical  layer  characteristics.  As 
a  result,  the  accuracy  (and  validity)  of  results  developed  under 
protocol  model  remains  unclear. 

As  discussed  in  Section  I,  physical  model  is  widely  accepted 
as  an  accurate  characterization  of  interference.  Unfortunately, 
although  physical  model  is  accurate,  there  is  much  difficulty  in 
carrying  out  analysis  with  physical  model  due  to  the  compu¬ 
tational  complexity  it  involves,  particularly  when  it  comes  to 
cross-layer  optimization  in  a  multi-hop  network  environment. 
As  a  result,  various  simplifications  have  been  employed  in 
recent  investigations.  In  [2],  Behzad  and  Rubin  studied  the 
special  case  that  the  same  power  level  are  used  at  each  node 
and  found  that  the  maximum  transmission  power  should  be 
used.  However,  for  the  general  case  where  each  node  can 
adjust  its  transmission  power  independently,  a  general  solution 
is  not  available.  In  [5],  Elbatt  and  Ephremides  proposed  a  two- 
step  approach  with  the  aim  of  using  a  minimum  power  vector 
while  supporting  as  many  users  as  possible;  routing  is  given 
a  priori  instead  of  being  part  of  the  optimization  problem.  In 


[3],  Chen  and  Lee  proposed  a  layered  (de-coupled)  approach 
for  QoS  scheduling  and  power  control;  while  in  [4],  Cruz  and 
Santhanam  proposed  a  two-step  approach  to  minimize  a  power 
cost  function  that  first  optimizes  link  scheduling  and  power 
control,  and  then  optimizes  routing.  Due  to  de-coupling  in  the 
solution  procedure,  these  approaches  only  yield  sub-optimal 
solutions. 

In  the  area  of  multi-hop  CR  networks,  although  there  has 
been  some  recent  efforts  (e.g.,  [17],  [18]),  none  of  them 
has  made  significant  progress  in  addressing  joint  optimization 
across  multiple  layers  via  the  physical  interference  model. 

On  another  line  of  research,  various  efforts  have  been  made 
to  study  asymptotic  behavior  (or  scaling  laws)  of  wireless 
networks  (see,  e.g.,  [8],  [12]).  These  efforts  differ  from  ours 
in  this  paper,  which  focuses  on  designing  optimal  cross-layer 
algorithms  for  finite  sized  network. 

VI.  Conclusion 

In  this  paper,  we  investigated  cross-layer  optimization  prob¬ 
lem  for  multi-hop  CR  networks,  with  joint  consideration  of 
solutions  at  physical,  link,  and  network  layers.  We  gave  a 
mathematical  characterization  for  power  control,  scheduling, 
and  routing  under  physical  interference  model.  We  developed 
a  centralized  solution  procedure  based  on  the  branch-and- 
bound  framework.  Using  numerical  results,  we  demonstrated 
the  efficacy  of  the  solution  procedure  and  offered  quantitative 
understanding  on  the  joint  optimization  at  different  layers. 
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